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\^ ' The two-body currents for s-wave pion photoproduction on the A^A^ system are derived to 0{Q'^) 

•^ ■ in chiral perturbation theory. For the interesting case of ^Si <-*• ^Sq transitions, we show that an ax- 

^^ ■ ial isovector two-nucleon contact term connects the short-distance physics of pion photoproduction 

to pion production and several important electroweak reactions. We also find that the standard 
\f^ ' chiral Lagrangian gives a ^irnNN vertex that have not been explicitly mentioned in previous lit- 
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^^ ■ erature. The corresponding Feynman rule is presented here and some processes where it should be 
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^D . important are briefly discussed. 
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Photo- and electroproduction of pions on nucleons have been intensely studied in the 
framework of chiral perturbation theory (xPT). Calculations exists for '-yd -^ ir^d [1], '~f*d —>■ 
Ti^d [2], 7*^*^A^ — > TiN [3], 7p -^ n'^'ii^p [4], and jd — ^ nnn^ [5]. Recently, also radiative 
pion absorption on the deuteron has been calculated using xPT [6, 7]. The latter work is 
motivated by the crucial importance the n^d -^ nwj measurements play for the extraction 
of the neutron-neutron scattering length (a„„) [8] and the corresponding need for a precise 
and accurate theoretical calculation. 

The advantage of using xPT (apart from providing a tractable low-energy limit of QCD) 
is that it is an effective field theory based on an expansion in the ratio x = Q/Mqcd, where 
Q ~ ?7i^ is the size of typical momenta and energies in the problem {m.,^ is the pion mass) 
and Mqcd = 0{1 GeV) is the energy scale at which this EFT breaks down. In the following, 
we will assume that Mqcd = M (the nucleon mass) and that the electron charge e < also 
counts as the small scale Q. (A detailed review of xPT can be found, e.g., in Ref. [9].) This 
expansion not only provides a prescription on how to estimate errors, but also suggest what 
needs to be done in order to improve the precision of the calculation. 

In this brief note we derive the 0{Q'^) two-body contributions to s-wave pion photopro- 
duction in Coulomb gauge. The amplitudes for neutral pion photoproduction were used in 
Refs [2] , though no expressions were given. We will find that the standard chiral Lagrangian 
gives an additional 'y-KiiNN vertex that have not been explicitly mentioned in previous lit- 
erature, though it was included in the calculations of [2, 10]. We derive the corresponding 
Feynman rule and discuss some reactions where it should be important. We will also ex- 
plicitly prove that the two-body contact term needed in [7] to reduce the theoretical error 
in the 7r~d ^ nw-f calculations is indeed the same as the one appearing in weak pp fusion, 
tritium (3 decay, the hep process (^Hep -^ ^Hee"'' Ue), i'{u)d breakup, and fi~d -^ nnun, as 
well as in p-wave pion production on two nucleons and the chiral three- nucleon force (3NF). 
This identity can be established only if the additional Feynman rule is included. 

The 7rN{N) terms of the chiral Lagrangian relevant for the 0{Q^) amplitudes are [9, 11] 

^il = N^tvD + gAS-u]N, (1) 
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where v is the nucleon four-velocity, D is a (chiral) covariant derivative, S is the Pauli- 



Lubanski spin vector, and gp^ parametrizes the unknown short-distance physics of the 
nucleon's axial-current matrix element. The low-energy constants (LECs) q = 0{-^) 
parametrize pion rescattering on one nucleon, while the di = Oi^jjp) parametrize the two- 
nucleon axial isovector contact term [19]. The latter have to be determined from fits to 
suitable two-nucleon data. Here u^ is an axial four-vector which contains the pion field, and 
when expanded takes the form: 
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where V^i (A^) is an external vector (axial) field and fj^ = 93 MeV is the pion decay constant. 
Note that only the isovector part of the vector field contributes at this order, (u" is obtained 
by factoring out r".) 

At next-to-leading order [NLO or 0{Q^)] a slew of diagrams appear, given in Fig. 1. 
These diagrams are the same as the ones in Refs. [1, 2], except that here only diagrams 
giving s-wave pions are shown. The diagrams have been grouped together according to their 
pion propagator structure. Thus, the first row [(a)-(d)] have a Coulomb-like propagator 
since the energy of the photon has to be transferred through the exchanged pion in order 
to produce a final, real pion, while row two [(e)-(g)] has an off-shell (zero energy transfer) 
pion propagator, and row three [(i)-(k)] has one propagator of each kind. The initial virtual 
pion in (e) and (i) is emitted in an s-wave (from £^^_,_^^^), which makes these terms 
proportional to the (in this case) very small energy of this pion. Thus these diagrams are 
pushed to higher order because of kinematics and will not be considered here. In addition, 
diagram (e) will vanish in the Coulomb gauge. The stretched pion exchange diagram (1) 
was included in [1, 2], where it was evaluated in time-ordered perturbation theory, but its 
evaluation depends on the way the attached wave functions are constructed, e.g., if they 
are energy-dependent or not. For this reason, we will not derive explicit expressions for this 
type of diagram. Finally, diagram (m) contains the two-body contact terms given by the 
LECs di and d2- Explicit expressions for the surviving diagrams are given in the Appendix. 

In principle, diagrams (c), (d), (f), (g), (h), and (k) should all have contact terms, i.e., 
they contain constant pieces that gets Fourier transformed into ^^^■'(r) in configuration space. 

3 

However, for '^5'i ^^ ^Sq transitions the (jg^s) contributions from (c) and (g) will partially 
cancel against (k), while the {^3) contributions from (d) and (h) will be completely canceled 
by a term from (f). This leaves the contact terms given by C3, C4 -|- -^ [from (f)], di, ^2 
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FIG. 1: Next-to-leading-order [©(Q^)] two-body diagrams for s-wave pion photo-production on 
two nucleons. Only one representative vertex ordering is given for each type of diagram. The black 



disc indicates an insertion from L 
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[from (m)], and the surviving g\ piece from (c), (g), and (k) as the only relevant ones. The 
contact terms linear in g^ will always appear (for ^Si ^^ ^Sq) with the combined coefficient 
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as explained in the appendix. The gj^ contact term appears only in the n d —>■ nn'y and 



7<i — > nnn^ reactions. The reduced (dimensionless) LECs q and di are defined through 

As claimed in [7] , the expression (4) is the same as the one appearing in pp fusion, tritium (3 
decay, and the hep reaction [11], neutrino-induced deuteron breakup [14], muon capture on 
deuterium [15], p-wave pion production on two nucleons, and the chiral 3NF [16]. Because of 
Fermi-Dirac statistics, the LECs di will always appear in the linear combination di + 2d2 [11]. 
The gauged C4 term [derived from the Lagrangian (2) and part of diagram (f)] has not, 
to the best of our knowledge, been explicitly mentioned in the literature before, i.e., it is 
missing in the Feynman rules derived in Ref. [9] [their Eq. (A. 31)] [10]. The Feynman rule 
for this term can be written (in a notation reminiscent of [9]) as 
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<5"V' - <5"V [5 ■ gi, 5 ■ e^] - U^'r' - 5"^" [5 • gs, 5 • e^] , (6) 



where a, 1(6, 2) annotate the incoming (outgoing) pion's isospin and four-momentum. The 
corresponding pion rescattering diagram (f) is the only way to obtain the C4 contribution 
needed for d in vr'rf -^ nn'^ or '-yd -^ nnn'^. In these ^Si ^-^ ^Sq transitions, the rescattered 
pion emitted by the gauged C3 and C4 terms is in an s-wave relative to the nucleon pair. This 
should be contrasted with 7*^*-'^ -^ iv'^d, where C4 (but not C3) can contribute and gives a p- 
wave pion only. This C4 term was included in the 'y*d — > ir^d calculations of Refs. [2, 10], but 
without being explicitly referred to. It will also contribute to 0{Q^) (charged) double-pion 
photoproduction on a single nucleon above threshold. 

Explicit expressions for the different surviving diagrams of Fig. 1 are given in the ap- 
pendix. As mentioned above, there appears an extra contact term from diagrams (c), (g), 
and (k) that did not appear in the other reactions mentioned. However, it enters with a 
known coefficient and thus does not destroy the linear correlation between the Gamow- Teller 
matrix element and the FSI peak height that was used in Ref. [7] to reduce the theoretical 
error in the extracted a„„. Even so, this term does shift the location of the correlation (it 
changes the FSI peak height) and is therefore necessary for an accurate extraction of a„„. 

In this paper we have derived the 0{Q^) two-body amplitudes for s-wave radiative pion 
capture or pion photoproduction on two nucleons. We have also shown that the axial 
isovector contact term for these processes can be parametrized by the same single constant 
that is used for certain electroweak reactions (like pp fusion and the hep process), p-wave 



pion production on two nucleons, and the chiral three-nucleon force. This was accomphshed 
through deriving an additional, previously not published, Feynman rule (parametrized by 
C4) for a photon coupling to two pions and one nucleon. This vertex should play a role in 
TX^d -^ nwy [7, 17] and 'yd -^ nrnr^, as well as in 7^*)^ -^ n^d [2] and 7A^ -^ txtxN beyond 
threshold, for the latter two reactions resulting inp-wave pions. Depending on selection rules 
for a particular reaction, recoil corrections might become important [5]. Further work on 
this may be needed — together with a calculation of C(Q^) single-nucleon amplitudes and an 
investigation of when the time-ordered diagram of Fig. 1(1) needs to be included. This will 
also necessitate a thorough understanding of how to consistently include 1/M corrections in 
operators and the NN potentials used to derive wave functions. We note that weak axial 
two-body currents and contact terms have been studied recently, from a somewhat different 
perspective, by Mosconi, Ricci, and Truhlik [18]. 

The results presented here are important for the calculation oi Ti~d —^ nn^ for two 
reasons. Firstly, it supports the conclusion of Ref. [7] that there is a correlation between 
the short- distance physics in electroweak reactions and n^d — > nw-f. Secondly, it gives 
expressions for the long-range 0{Q'^) currents that should give a correction to a„n- These 
currents, together with consideration of the issues raised in Refs. [5, 18], will be included in 
forthcoming work on a precision computation of the 7r~d ^ nw-f spectrum [17]. 

It is a pleasure to acknowledge discussions with Fred Myhrer, Kuniharu Kubodera, and 
Daniel Phillips. I appreciate communications from Veronique Bernard and Ulf-G. MeiBner 
regarding their use of the Feynman rule presented here. This work was supported in part 
by the National Science Foundation grant PHY-0457014. 

APPENDIX: EXPRESSIONS FOR O(g^) TWO-BODY DIAGRAMS 

We give here momentum-space expressions for the diagrams of Fig. 1 for the case of 
s-wave pion photoproduction. These formulas can also be used for 7r~d ^ nw-f by time- 
reversal invariance. In the equations below, we assume that p = p" — p' and 2P = p" + p', 
where p'(p") is the relative momentum of the incoming (outgoing) nucleon pair. Also, e^ 
and e^ are the deuteron and photon polarization vectors, while k is the photon wave vector 
and /to(i) is the isoscalar (isovector) anomalous nucleon moment. Everything is expressed 
in the overall rest frame. When doing the replacement (1 ^^ 2), the momenta p and P 
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change sign. It is assumed that contact term contributions should not be included when 
evaluating the expressions (A.l)-(A.ll) below (i.e., that ^^'^•'(r) terms should be dropped 
when performing the Fourier transform). They are instead given explicitly in Eq. (A. 12). 

1. Coulombic diagrams 

The momentum-space expressions for the first row of diagrams in Fig. 1 are 
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2. Non-Coulombic diagrams 



The momentum-space expressions for the second row of diagrams in Fig. 1 are 
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3. Double-propagator diagrams 



The moinentuin-space expressions for the third row of diagrams in Fig. 1 are 
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4. Contact term 



The contact term that contributes to ^Si ^-> "'^S'o transitions at 0{Q'^) is given by 
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Obviously, this contact term can not contribute to processes involving a neutral pion. When 
put into a ^Sq —> ^Si matrix element, specialized to -K'd — > nw-f (or rather 'jnn -^ rfvr"), 
and summed over nucleon spins and isospins, this operator reduces to 
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where d is given by Eq. (4). The g\ term is the extra piece coming from diagrams (c), (g) 
and (k). 
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